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We study a modified Bullard dynamo and show that this system is equivalent to a nonlinear 
\£i ' oscillator subject to a multiplicative noise. The stability analysis of this oscillator is performed. Two 

bifurcations are identified, first towards an "intermittent" state where the absorbing (non-dynamo) 
state is no more stable but the most probable value of the amplitude of the oscillator is still zero and 
£N) ' secondly towards a "turbulent" (dynamo) state where it is possible to define unambiguously a (non- 

zero) most probable value around which the amplitude of the oscillator fluctuates. The bifurcation 
diagram of this system exhibits three regions which are analytically characterised. 

£N| ■ Introduction 



In this paper, we aim to study how a bifurcation is modified when the control parameter fluctuates. This is an 
important question regarding the problem of self-generation of a magnetic field by means of a conducting fluid, known 
as the dynamo instability. From the experimental point of view, a working dynamo is one where the fluid is turbulent 
and thus the velocity (which is the control parameter of the dynamo instability) fluctuates. However, experiments 
have been designed using the mean-field (mainly time averaged velocity field) and the effect of the fluctuating part 
of the velocity field still remains an issue Q, Q ■ The main framework to study the effect of a fluctuating velocity 
field on the dynamo threshold is that of "Mean-Field theory" (MFT) of dynamo 0, Q where the induction equation 
£NJ ■ is supplemented by two terms, an "a effect" which corresponds to the creation of large-scale magnetic field by helical 
\ small-scale motions and a u j3 effect" which accounts for a turbulent diffusivity. The validity of the mean- field theory 
has however been questioned due to the statistical nature of the turbulent dynamo theory: due to fluctuations, there 
may be some additional effects which MFT cannot account for [3, Ej ■ Another limitation in the MFT approach is due 
to the fact that it is purely linear, disregarding the effect of the magnetic field on the velocity field, and consequently, 
may overlook some effect due to the interplay between noise and non-linearity, such as stabilisation by noise 0, E| , 
ly-j , noise-induced transitions^^]], stochastic resonance ^ij, etc. 

Based on the fact that the fluctuating parameter (the velocity) multiplies the magnetic field, it has been argued 
""^3 \ that the magnetic field above threshold should be intermittent, exhibiting sporadic growth of the magnetic field 
• i-H ' [l2L H^ l . This feature has been evidenced in the case of a stochastic model of the turbulent dynamo 0] but only 
under simplified assumptions: zero magnetic viscosity and 5-correlated noise. Consequently, this is not clear if a "real 
• • dynamo" would exhibit intermittency. To shed some light on this problem, we study the much simpler problem of a 
solid dynamo with a scalar control parameter and we show that this problem can be mapped into that of a nonlinear 
oscillator with a fluctuating frequency, a problem that has been widely investigated in the literature. The question 
of its stability has been questioned both in the linear case 0] and the nonlinear one 0, llCjL IT7L fl8| . We then use 
known results to draw the stability diagram of this system and show that, in the non-linear case, the system bifurcates 
towards an intermittent state (characterised by a mostly null field, except for some small amount of time where the 
field exhibits bursts of activity). We also show that if the intensity of the fluctuations is increased, that particular 
behaviour can disappear and the field oscillates around a well defined mean value. 

After introducing the disk-dynamo model, we derive the analogy with a noisy oscillator (section [HJ) and in section 
IIIII both the frontier between the trivial state and the intermittent state (already calculated by 0|) an d between the 
intermittent state and the oscillatory state are analytically calculated. We also address the effect of noise correlation 
by means of numerical simulations (section HVfl . 



I. MODEL 

The Bullard (or homopolar) dynamo [nil ] is the first example of a magnetic instability triggered by a conductor in 
motion. The experimental device is depicted on figure ^ a conducting disk rotates around its axis at the angular 
velocity u> and a small magnetic field is applied in the vertical direction inducing a current in the disk from the axis 
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towards the edge of the conducting disk. Then this current flows in a conducting loop which, if orientated in an 
appropriate way, induced a magnetic field orientated in the same direction as the initial one. Thus an infinitesimal 
magnetic field can be amplified by this mechanism, leading to a dynamo. Here we study a modified version of the 




FIG. 1: The homopolar (or Bullard) disk-dynamo. 

original Bullard dynamo where an azimuthal current at the edge of the disk and mechanical friction are permitted. 
The azimuthal current has been introduced by Moffatt 20] who showed that the hypothesis of a purely radial current 
would violate the conservation of magnetic flux in the case of a perfectly conducting disk. According to Hide [2l| . 
neglecting the mechanical friction is unwarranted for it would render the disk-dynamo structurally unstable. Three 
equations are enough to describe this system [22, Hj| : 

x = q(y~ x) 

y = xZ + mx — (m + l)y (1) 
Z = g[l — (m + l)xy + mx 2 ] — fZ 

where Z is the dimensionless angular velocity of the disk and x and y are the two components of the magnetic flux, 
across the disk and across the loop (see [2^ for the precise meaning of the different constants). The two first equations 
represent the effect of the rotating disk on the magnetic field and the third equation is the evolution of the angular 
velocity, subject to Joule heating and externally applied torque (the term proportional to g) and a mechanical friction 
(the term proportional to /). A linear instability analysis around the solution Z — Zq (given angular velocity) and 
x = y = (zero magnetic flux), proves that the system becomes unstable as soon as Zq > 1 . 

II. A NON LINEAR NOISY OSCILLATOR 

To mimic the turbulent dynamo problem evoked in the introduction, we make a kind of kinematic approximation: 
we prescribe a fluctuating velocity field and study the generation of the magnetic field induced by this given velocity. 
Practically, this accounts to disregard the third equation and study the two first equations only. However we intend 
to compare the result of this procedure to the numerical analysis of the full system, and while doing such, it could be 
important to keep track of the non linear term in the third equation (the back reaction of x and y on the velocity Z) . 

Specifically, we decompose the velocity field as a sum of a mean part and a fluctuating part: Z = Zq + T(t). For 
simplicity, we assume that the fluctuating part of the velocity field is white in time : (T(t)r(t')) = 2D5(t—t'). Applying 
this procedure to the two first equation of (Q, we are left with a linear stochastic system. It is well known that it is 
difficult to define unambiguously a threshold for this system and it is necessary to account for the non-linearity induced 
by the third equation. In consequence, we introduce a non-linear term and write Z — Zg+T(t) + [— (m+l)xy+mx 2 ]g / f 
to mimic the feedback of the variables x and y on the intensity of the dimensionless angular velocity. 

The constitutive equation of this system can be transformed into the equation describing a nonlinear oscillator in 
presence of noise and damping: 



x + (1 + (3\x 2 )x -ax = - Ax 3 , with (£(*)£(*')) = A< ^ - f ') . 



(2) 
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where the time and the variable x have been adimensionalised by (to. + 1 + q) . The parameters are: 

q(Z - 1) A 2Dq 



a 



A = 



(m + l + g) 2 

gg/f 

(m + 1 + q) A 



A = 
and 



(to + 1 + qf ' 

_ (TO+l + g )(m + l) 



(3) 



The equation J2J) is similar to the equation of a Duffing oscillator except for the presence of non-linear terms. The 
undamped version of this oscillator has first been studied in the context of Anderson localisation problem [24| and 
more recently introducing a damping and studying the linear stability of this oscillator [Isj . Here we focus on the non- 
linear version of the oscillator and show that the transition observed in the non-linear system is of a very particular 
nature and that it is possible to identify a series of two bifurcations in this system. To study @ , we make the following 
change of variables from the Cartesian coordinates (x,x) to the polar ones (r, 9): 



(4) 
(5) 



x = r cos(#) and x — r sin(#) . 
Then, the evolution equations for the new variables r and 9 can be derived: 

f = r sin 9(cos9 — sin 9) + ar sin 9 cos 9 + r cos 9 sin 9£(t) 

- Ar 3 (cos 3 9 sin 9 + cos 2 9 sin 2 9) 
9 = - sin 6>(cos 9 + sin 9) + a cos 2 9 + cos 2 9£(t) 

-Ar 2 (cos 4 6> + /3cos 3 <9sm60 . 

We are left with a system of two stochastic equations for r and 9. From these equations, using standard techniques 
of stochastic process [25|, it is easy to write a Fokker-Planck equation which governs the evolution of the density 
probability H (t, r, 9) to observe a couple (r, 9) at time t: 



d t H = C r H + C e H + KH , 



(6) 



where C r and Cg are two differential operators involving respectively derivative with respect to the radial variable 
and the angular variable: 



- [a sin 9 + sin 6»(cos 9 - sin 9) + — cos 2 6»(cos 2 9 - sin 2 9) d r (rH) 
A(cos 3 9 sin 9 + f3 cos 2 9 sin 2 9)d r (r 3 H) + ^ cos 2 9 sin 2 9d r [rd r (rH )] , 



(7) 



do 

A 
~2 



(sin6*(cos0 + sin 6) -a cos 2 9 + Ar 2 (cos 4 9 + (3 cos 3 9 sin 9)) H 
d 9 [cos 2 9d e (cos 2 9H)} , 



KH = Ad r dg[rsm9cos 3 9H] 



A. Radial density of probability 

The partial differential equation 10 admits simple analytical solutions with separation of radial and angular vari- 
ables: H(r,9) — P(r)G(9). Then, if we integrate equation 10 relatively to the angular variable, we obtain the 
following equation for the probability density of the radial variable: 

d t P = ad r [rd r (rP)] - bd r [rP] + cd r [r 3 P] , (8) 

with the following parameters expressed through an average over the angular variable ((•)# = / • G(9)d9): 

a = j (cos 2 9 sin 2 9) e , (9) 

b = (a sin # cos 6* + sin6'(cos6 | — sin 6*) + — cos 2 #(cos 2 9 — sin 2 9))g , 
c = A(cos 3 6'sin6'-l-/3cos 2 6'sin 2 6') e . 
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It is then easy to find a stationary solution of equation JSJ: 

p.( r ) = I r »/a-l exp[ _^] with Z =\ 
Zj la I 

When b is negative, one can check that the distribution is not integrable in zero. In that case, the only admissible 
solution is a Dirac function centred around zero: P s (r) = S(r) which corresponds to a solution x = y = at long 
time. Indeed, one can check that it is always a solution of equation J2J). Then, we can identify two bifurcations 
whether the control parameter is taken to be the mean value of r or its most probable value. This scenario as already 
been evidenced in the case of a stochastic modelling of the dynamo effect [l^. To compute the threshold value 
corresponding to this two bifurcations, one needs to characterise completely the probability density l|10|) and thus to 
calculate the coefficient defined above. To achieve that, one needs to compute the stationary probability density of 
the angular variable. 



. r. 



(10) 



B. Angular distribution of probability 



To express the coefficient a, b and c, we have to compute the distribution probability of the variable 0. This can 
be done by averaging equation JHJ over the r variable. This leads to the following equation: 



d t G = -d e [(acos 2 0-sin0(cos0 + sin0))G] + y<9 e [cos 2 0<9 e (cos 2 0G)] 
+ ma e [(cos 4 + /?cos 3 0sin0)G] , 



(11) 



where we set m = A J r 2 P s (r)dr = bX/c. It is interesting to notice that this coefficient does not depend on A and thus 
that the probability distribution does not depend on the intensity of the non-linear terms. It is convenient to solve 
the preceding equation to introduce the variable z — tanO — x/x which leads to the following Fokker-Planck equation 
for the probability distribution I(z,t): 



dj = -a 



l + (3z 



.1 + z 2 l + z 2 {1 + z 2 ) 
This equation can be solved in the stationary case {ptl — 0) to obtain: 



2 z 



I(z) 

with 



1 

Z 
at 



exp 



— oo 
2 



3 



dz 



(12) 



(13) 



arctan(t) + ^log(l + t 2 ) 



This expression is identical to that of |18| except for the last term which comes from the nonlinearity and which is 
important for the structure of the oscillator bifurcation diagram, as we now discuss. 



III. BIFURCATION SCENARIO 



The bifurcation of the noisy oscillator is of a rather complex nature. Indeed, depending on the control parameter, 
two thresholds can be identified. This is to be contrasted with the case without noise where a < corresponds to an 
absorbing state where the oscillator relax towards x = and a > where the solution at long time is x — -Ja/X. We 
now characterise more deeply the nature of the bifurcation when the noise intensity is not zero. 



A. b > 0: bifurcation towards an intermittent state 



When b becomes positive, the distribution (|10[1 is integrable near the origin and can thus represent a meaningful 
probability distribution function. However, the distribution diverges in zero, an indication that the most probable 
value is still zero. This is characteristic of an intermittent state where the signal exhibit bursts of activity separated by 
quiescent epochs where the norm is close to zero. To illustrate such a behaviour, we performed numerical simulations 
of the equation J2J and a typical snapshot is shown on figure [3 
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FIG. 2: Result of a numerical simulation of system J5J with A = 2, A = 1, /3 = 3. The left panel is for a — —0.2, corresponding 
to an absorbing state and the right panel for a = 0.2, corresponding to an intermittent state. 



On the left hand side, we see that the variable x after a short transient time (compare the time interval of the two 
snapshots) relax towards the absorbing state. Increasing the parameter a (keeping the other parameter fixed at the 
value of figureEl), one notice that a first bifurcation occurs (around a = 0.17) which leads to a state where the variable 
x is most of the time close to zero but exhibits bursts (cf the right hand side of the figure). We call this type of 
behaviour intermittency. 

To characterise the intermittent state, we use the time series of x to compute the different parameters that appeared 
in the previous section a, 6, c and m. Once this parameters are known, we can compute the probability distribution of 
the variable r and z given by the expression l|10|) and (|13fl . The theoretical distribution for the angular variable shows 
a good agreement with the numerical ones as is shown on figure Concerning the radial distribution, the agreement 
is good for low values of r but not for high values. This is because the decoupling approximation is supposed to be 
relevant for low energies (or low r). However, because we are interested in the stability of the r = fixed point, the 
approximation is good in the region of interest. 




FIG. 3: Comparison of the simulation with A = 2, A = 1, /3 = 3 and a — 0.2 corresponding to an intermittent sate and the 
analytical prediction (in straight line) for the probability density of the variable z (left hand side) and r (right hand side in 
log-log scale). The parameters have been found to be numerically a — 0.14, b = 0.019, c = 0.47 and m = 0.04. 



B. b > a: bifurcation towards a turbulent state 



When b becomes larger than a, one can easily check that the distribution has a well defined maximum for r = 
\J '{b — a) I 'c. A typical example of this behaviour is shown on the left of figure 0] 
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FIG. 4: Time series (left) of the variable x and probability distribution (right) of the variable r for values of the parameters 
corresponding to a turbulent state: A = 2, A = 1, /3 = 3 and a = 2. The numerically calculated value are: a = 0.088, b = 0.54 
and c = 0.25. 



In this new regime, the variable x(t) fluctuates around a well defined most probable value which can be seen both 
on the time series (left) and on the probability distribution which has a well defined maximum. One sees that the 
theoretical distribution does not fit the numerical one accurately anymore, a possible indication of the breaking down 
of the separation hypothesis. This prevents us to make any prediction on the mean value of a: in the turbulent regime. 
However, due to the fact that the theoretical distribution fits well the numerical one in the intermittent regime, it 
is still possible to calculate the boundary between the intermittent and the turbulent regime and to draw the phase 
diagram of the system. We will now proceed to this determination. 



C. Stability Diagram 

In previous sections, we identified a scenario of two bifurcations for the non linear oscillator: starting from the 
absorbing state x = 0, a first bifurcation leads to an intermittent regime (see figurc[3J) followed by a second bifurcation 
leading to a turbulent regime. These two regimes can be identified via the shape of the probability density of the 
variable r, or, equivalently, by means of conditions on the coefficient a and b. We now write these conditions in terms 
of the initial variable a, A, A and 0, so as to fully characterise the phase space. 

For the first bifurcation (towards the intermittent state), we notice that the position of the bifurcation line does 
not depend on the non-linear terms. Indeed, under this threshold, the only solution is x = and thus the non linear 
term is negligible. This has already been stressed by |l8j and we here just outline the main lines of the computation. 
First, multiplying equation {TJ by log(l + z 2 ) and averaging over z, one finds the following relation: 



z + Bz 2 . 

{z + m {i + z 2 y h (14) 



which reduces to b = (z) upon discarding the non linear terms. This is the expression for the Lyapunov exponent 
found by p4j in the case of Anderson's localisation problem and later by jTsj when studying the linear version of 
our oscillator. Using, the expression for the probability distribution of z H13fl . one may then express the condition of 
instability b > and find the bifurcation line a* (A). 

The condition for the appearance of intermittency can thus be written as b = (z) > 0. As shown by |Ts[ . this can 
be written: 

f +OC j r- \ 2 t< 1 , u3 \] f + °° du \2f. 1 w3\i 
J du ^ ex nA\ {a+ 4 }U -T2\\ > J ^ eXP [A{ {a+ l )u -u)\ ■ (15) 

The numerical integration of this condition give the bifurcation line a(A) which is drawn on figure|S]in full line. This 
curve delineates the parameter space between the absorbing state (or the "no dynamo" state), corresponding to x 
going to zero for long time and the (intermittent and turbulent) dynamo states. 

To draw the frontier between these two states, one need to re-express the condition b > a in terms of the initial 
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parameters. This can be written: 



/ z+0z 2 

z + (3z 2 \ A / z 2 \ 6A \ Z + W (7+?F / 
z + ™7^ 9\9/ > TT\7i 9V9 / wnere W = — = 7 -, (16) 



(1 + z 2 ) 2 / 2\(l + z 2 ) 2 / c / z±§£_\ 



where the brackets correspond to a mean using the distribution 1)13(1 . Notice that this condition and the distribution 
only depend on a, A, j3 and m. Furthermore, it is easy to see that the factor A cancels out in the expression for m. 
As a consequence, the bifurcation line a(A,/3) between the intermittent and the turbulent state does not depend on 
the particular value of A, i.e. on the absolute intensity of the non linear terms. This is a rather remarkable feature of 
the bifurcation. The result of the numerical integration of the two equations Ijlfcifl is shown on the figure |SJ One sees 
that the evolution of this threshold is monotonous: when the noise is increased, the transition from the intermittent 
to the turbulent state is delayed. On the contrary, the transition from the absorbing state to the intermittent one 
is first increased for weak intensities of the noise (stabilisation by noise) and then is lowered below its deterministic 
value for more powerful noise (a reentrant transition that has been pointed out by |l£j|). 




IV. INFLUENCE OF THE NOISE CORRELATION 



The next step will be to compare the stability diagram of figure with a numerical simulation of the system l|T]l. 
Indeed, the white-noise hypothesis has been introduced for the sake of simplicity but it is really far from clear if a 
turbulent velocity has something to do with a white noise. In [l^, the effect of an exponential correlation of the 
noise (Ornstein-Uhlenbeck process) has been studied. It has been showed that the system behaviour was qualitatively 
the same: a stabilisation by the noise for weak noise intensities and a noise induced bifurcation for stronger noise. 
However, we know also that a turbulent velocity has not a pure exponential correlation function. 

In order to get an insight into the influence of a "real noise", we will use data from a fluid dynamics experiment. 
We will here use a signal of velocity from a von Karman experiment. Details about the experimental setup can be 
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found in |26( but for our purpose, it is enough to know that this device gives rise to a fully developed turbulence. On 
figure IHI we show the correlation function that can be extracted from the fluctuating part of the velocity signal: 

i r T 

c(t) = — v' \u)v' (u + T)du T->oo. (17) 

J- .In 




FIG. 6: Correlation function of the turbulent velocity from a von Karman experiment. 

From the correlation function, we can define the parameters that will serve us to compare the correlated case, to 
the calculation we have made in the previous section. We will use the following definitions for the intensity of the 
fluctuation and the correlation time of the signal: 



A = 



C(t)dt 



and 



C(0) 



(18) 



From the turbulent signal, we found a correlation time of r ~ 0.14 s. By multiplying the turbulent fluctuating velocity 
v'(t) by a suitable quantity, we can obtain a signal with any value for the parameter A, while keeping the correlation 
time constant. Then, we use this signal in equation J5J as a noisy term = x' (t)). Using the experimental noise, 
we observe as previously a succession of two bifurcations, first towards an intermittent state and then to a fluctuating 
one. However, the second transition is much more smooth than in the (5-corrclatcd case and is thus very difficult to 
identify. We will not discuss more about this transition. On the contrary, the first transition, from the absorbing state 
to the intermittent, is still easy to identify. We thus perform some numerical simulation, keeping A at a constant 
value and increasing the control parameter a. The bifurcation line, separating an absorbing state (under the curve) 
and an intermittent state (above) is compared to the (5-correlated case on figure 

We see that the qualitative behaviour (stabilisation for weak noise and destabilisation for higher intensities) is not 
altered by the correlation of the noise. Moreover, we see that this scenario is reinforced: both the stabilisation and 
the destabilisation of the instability is much more marked than in the previous case. 



SUMMARY AND CONCLUSIONS 



We studied a modified Bullard dynamo model showing that it can be mapped onto the problem of the stability of 
a non-linear oscillator in presence of a 5-correlated noise. We studied the probability density of the position of the 
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FIG. 7: Stability curve for the 5 correlated model (straight line) and for the "experimental noise" (stars) . 



oscillator induced by the noise and showed that a non trivial scenario of bifurcation emerges because of the noise. 
First, looking at the mean value of the position, we found a transition from an absorbing state where the oscillator 
is trapped at his rest position for long time, to an intermittent state, where it stays close to the absorbing position 
for long time but sometimes makes excursion far from this position. However, the most probable position for the 
oscillator is still at the rest position. Thus, using the most probable value as an order parameter, we showed that a 
second bifurcation takes place for higher intensities of the control parameter. The evolution of the two bifurcation 
threshold has been studied relatively to the intensity of the noise. A recent result by [lg has been recovered: for 
weak intensities, the instability is delayed whereas for higher noise, there is a reentrant transition induced by the 
noise. Furthermore, we showed that it was not the case for the threshold of the second transition which increases 
monotonously with the strength of the noise. Then, we checked what could be the influence of the correlation of 
the noise on this bifurcation scenario and showed that the qualitative picture was not changed but the threshold 
displacement was amplified by these correlations. 

The transition to an intermittent state, known as the "on-off intermittency "in nonlinear physics has regained 
some attention in the past year because, despite its generality in chaotic system and stochastic system driven by 
multiplicative noise, there is very few experimental evidences of such a behaviour. In a recent paper, Aumaitre et at 
|27j have shown that the existence of intermittency was mainly monitored by the low-frequency components of the 
noise and because of filtering, these components may not be present in experimental devices. Regarding the dynamo 
effect, it is still a question if the first bifurcation will be towards an intermittent state ? The preceding study show that 
the intermittent bifurcation observed by [14;] is not an artefact's due to the white noise approximation. However, this 
behaviour could still be induced by other approximations such as the zero diffusivity limit. Furthermore, in a realistic 
dynamo experiment, there is other effects that could prevent us to observe a bifurcation towards an intermittent state, 
e.g. the presence of the external Earth magnetic field or the back-reaction of the magnetic field on the flow. 

To study the effect of this back-reaction, we started the numerical study of a "dynamical" model, i.e. where the 
control parameter is not prescribed statistically but evolves with a chaotic equation of motion. In the kinematic 
approximation (where the effect of the magnetic field on the velocity field is neglected), the evolution of the magnetic 
field will be governed by the two first equation of with a chaotic variable Z. We aim to compare the results from 
the numerical study to the model presented here and investigate, firstly if the back-reaction of the magnetic field 
change anything in the bifurcation scenario and secondly what is the effect of a real correlated noise compared to the 
white one. This work is under progress. 
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